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\ We derive the so-called Barbieri-Remiddi solution of the Bethe-Salpeter equation 

in QED in its general form and discuss its application to the bound state energy 
spectrum. 
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1 Introduction 



X 

The Bethe-Salpeter (BS) equation is usually considered the rigorous framework in which 
to approach the bound state problem in QED. The increasing precision of the experimental 
data concerning QED bound states (e. g. decay rate, energy levels, etc. for some recent 
reviews we refer the reader to 0) makes more and more urgent to effort the evaluation of 
physical quantities by handling the BS equation with a systematic and unified formalism. 

In this paper we will focus our attention to the bound state energy levels in QED and will 
discuss the so-called Barbieri-Remiddi (BR) formalism. This formalism was first suggested 

for positronium H £§] (for a clarifying quanto-mechanical example see also ||), but has been 
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used in recent years also for hydrogenic atoms 0, QCD bound states ]7) and scalar-scalar 
bound states M . The main idea is to solve exactly the BS equation for a suitable zeroth-order 
kernel containing the relevant binding interaction (i. e. the Coulomb potential) and then to 
perform a perturbative expansion in terms of the difference between the complete two-body 
kernel and the zeroth-order one. What is appealing in this approach is that the zeroth-order 
solution is completely known in analytic closed form. Therefore the perturbative expansion 
obtained in this way is completely self-contained and does not need to be improved for higher 
correction in the fine structure constant a. 

In the following the BR formalism will be derived in the general case of muonium (i. e. 
different masses). This result is new and contains the positronium and hydrogenic case as a 
particular one. Moreover it furnishes a way to treat radiative and recoil corrections in the 
same theoretical framework, which seems to be very promising. 

The paper contains two main sections. In section 2 we derive the perturbative expansion 
of the energy levels from the BS equation in the so-called Kato formalism. In section 3 we 
derive in some detail the BR solution for muonium. Section 4 is devoted to some comments 
and conclusions. 

2 The Bethe— Salpeter equation 

In this section we review some basics concerning the Bethe-Salpeter equation in QED and set 
up the theoretical background for the next section. The main result will be the perturbative 
expansion of the energy levels of the two fermion bound state given by Eq. ([12]). 

Let us consider a system of two fermions (of mass mi and m 2 and electric charge — e and 
Ze respectively) like muonium. The four point Green function G is the sum of the Feynman 
graphs shown in Fig. [I] (notice that for a particle-antiparticle system, like positronium, one 
has to add the annihilation graphs). Let us define G the two fermions free propagator: 

G (E;p) = s£\E 1 ;p)S i *\-E 2] p) T J S F (E j;P ) = ' — , (1) 

Eaj v + p — m + ie 
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Figure 1: Four point Green function G: graphs up to two loops. 



where ^ and ^ refer to the two fermion lines, Ej = E rrij/ (mi + TO2), E is the bound state 
energy and p, q are the relative momenta of the outcoming and incoming particles in the 
centre-of-mass reference frame []. It is found that the Green function G satisfies the equation: 

d A k 



G(E;p,q) = G (E;p) 



(27r) 4 5 4 (p- g ) + 



(2*) 



-K(E;p,k)G(E;k,q) 



(2) 



This equation (for simplicity we have neglected the spinor indices) is known as the Bethe- 
Salpeter equation ffl. The kernel K, describing the interaction between the two fermions, 



1 Let p^ and p( 2 ' be the momenta of the outcoming particles, and 

rj (11 , (2) _ A* (1) M (2) 

P — p y > + p y ' , p — p y ' p y ' 



mi 



rri2 



with /i = mimz/ (mi +TO2) the reduced mass of the two particles. In the centre-of-mass frame p W = — p ' 2 
implies P = (p^ +Pq 2 \o ) = (E, ). From the previous equations we obtain 



P 



(i) 



mi 



V 



(2) 



-P 



mi + m 2 
m 2 



P=( Po + Ei,p), 
P = (-po + E 2 ,-p), 



mi + m 2 

where Ej — E mj /(mi + m 2 ) and therefore 

E = Ei + £-2- 

Finally, we note that in the static limit 

E — > mi + m 2 , Ej — > 
In the same way we can treat the incoming momenta. 



is not known in analytic closed form and is given by all the two particle irreducible graphs 
without external legs shown up to two loop in Fig. |2|. Graphically the BS equation can be 
represented as in Fig. |^. 



K 




Figure 2: The interaction kernel K . 



G, as a function of E, has simple poles in the bound state energy levels E n _ || (n... 
is a convenient set of quantum numbers which classifies the levels). Since the Coulomb 
interaction is contained in K, E n _ (without mass terms) has to coincide with the Bohr 
levels at the leading order in a. Therefore for any n..., we can write 

G(E) = + G n ..{E) , E n ... = mi + m 2 - £ (— )' + ■■■ , (3) 

E — E n ,„ I \ n J 

where R n ... is the residuum at the pole, G n ... is non singular in the limit E —>■ E n ___ and \x 

is the reduced mass of the two particle system. From now on we will neglect the explicit 

indication of the momenta in the argument of the functions, where considered not strictly 

necessary. 

Inserting (|3|) into the BS equation and comparing the residua, we obtain: 

R n ... = G (E n ...)K(E n ...)R n ... } (4) 

which is known as the homogeneous Bethe-Salpeter equation. Moreover, from the compari- 
son of the non singular parts in E — E n „ we obtain the normalization condition [ ID |: 

R n ... = Rn... (Go 1 - K) (E n J R n ... . (5) 



Figure 3: Inhomogeneous BS equation: G(E) = G (E) + G (E)K(E)G{E). 



The BS equation (|2|) is, up to now, not solvable in analytic closed form. Let K c be an 
interaction kernel satisfying the following two properties: 

i) K c reproduces the same non relativistic limit of K, i.e. the Coulomb potential V = 
—47rZa/p 2 times some spinorial factors; 

ii) the BS equation for K c : 

G C {E) = G (E) + G (E)K C (E)G C (E), (6) 
is analytically solvable in closed form. 

With these assumptions it is possible to solve the BS equation for G at least perturbatively 
in terms of 5K = K — K c and to give a perturbative expansion for the bound state energy 
levels (the poles of G). 

From the property i) it follows that G c has simple poles for E = E^. These poles are 
surely more degenerate than the poles of the complete Green function G and give back, at 
the leading order in a, the Bohr levels: 

G C (E) = + d c n (E) , £n = ™i + ™ 2 ~(v) 2 + (7) 

where R c n is the residuum at the pole and G c n is non singular in the limit E —>■ E c n . The 
sum is extended over all the degenerate states for each nsN. The residuum R c n satisfies 
the analogous of equations ([|) and (||): 

J£„ = Go(E c n )K c (E c n )Rl__ , (8) 

K... = K... §e (Co 1 - O (K) K... ■ (9) 
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From the definition of SK and from (|2]) and ([]) we obtain the perturbative expansion of 
G in terms of 5K: 

G(E) = G C (E) + G C (E)SK(E)G C (E) + G C (E)5K(E)G C (E)5K(E)G C (E) + ■■■. (10) 

In order to obtain from (|10D the perturbative expansion of the poles E n „. we will use the 
so-called Kato perturbation theory JO]. Since the energy levels are the poles of G we can 
write: 

<f dz z Tr {G{z)O n } 

En... = , (11) 

i dzTi{G(z)O n } 

where T n is a closed curve in (E which contains only the poles E n _ and E° of G and G c 
respectively, O n _ is an operator which does not vanish on R n and R c n and Tr means the 
trace over the spinor indices. A convenient choice is 

On... = ^ (G - 1 - K c ) {ED it£... ^ (G- 1 - JT C ) 

Inserting (|T0|) in (0) integrating in z and taking into account we obtain (up to order 
«f 4 ): 

>'<■- >'„... — 

+ -L (5K(E^G c (E c n )5K(E c n )G c (E c n )5K(E c n : 

+ -d\- (mK)Gc(E c n )5K(E c n )) n ^5K(K 
d 



£>L. \dE 

+ 7^(^n))n.. 
n... 



I n I 



where 



+ lj^m(K)) n J 2 (^K(K)) + 0{5K% (12) 



D n _ is the degeneracy of the E = E n ,„ level and is defined to be 

D n ..=^(G^-K c )(E^ . (13) 

Equation flT!2| ) expresses the bound state energy E n .„ as an expansion in 5K. Since 8K is 
the difference between the sum of the infinit series of Feynman graphs drawn in Fig. |2| and 
the kernel K c , 5K is not known in closed form. Each graph of Fig. || contributes to ( |l2l) with 
a series of powers of a, because the dependence of the residuum R c n on the fine structure 
constant (like in the well-known non-relativistic case where the hydrogen wave-functions 
depend on a). For consistency with i) the explicit calculation must exhibit that to an 
increasing order in 5K it corresponds an increasing leading order in a in the contributions to 
the energy levels. In this sense expansion (|T2|) can be interpreted as a perturbative expansion 
in the fine structure constant. Once K c is explicitly given and the corresponding BS equation 
is solved (this means we have an analytic expression for E^, R c n and G^) the expansion ( p!2|) 



allows to obtain without any ambiguity the energy levels up to a given order in a for all the 
two-fermions bound states in QED. We emphasize that, in absence of an exact solution of 
Eq. (H), expression (|12|) could be evaluated only for an approximate choice of G c to improve 
at any increasing of the requested precision. 

3 The Barbieri— Remiddi solution 

In this section we work out with some detail the so-called Barbieri-Remiddi solution of 



equation @ (for and exhaustive description see [|12])- With this name we mean a zeroth 
order kernel K c which satisfies the previous given properties i) and ii) (in other words K c 
should describe correctly at the leading order in a the bound state and make solvable the 
corresponding BS equation (||)) as well as the solution of the corresponding BS equation. 
The BR solution was first given for the positronium 0. In the following we will give the 
generalization of that solution for a bound state of two fermions with different masses, i. e. 
muonium. Once K c is given, we solve the equation for G c and work out the poles E^ and 
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the residua . At that point the perturbative expansion of the energy levels ([L2|) will be 
completely defined. 

Let us define the energy projectors 

A±(p, rrij) = ; (14) 



with E p j = \jp 2 + rrij. In terms of A± the free fermion propagator Sp can be written as 

o ( p . „x / A+(p,m i ) 7 ° A_(p, m j ) 7 ° ^ 

b F {Ej,p)-% -5—= „ — + 



s p° + - £ pi + ie p° + Ej + £ pi - ie / 
Moreover, we define 



\(- -) = ( 16E pl E p 2E ql Eq2 \ 2 , , 

^ - \{E pl +m l ){E p2 + m 2 )(E ql +m 1 ){E q2 + m 2 )) 1 ' 

( 1+7° \ (1) / I-7 \ (2)T 

x ( 7°A+(p,mi)— - — A+(g,mi)j (7°A_(g,m 2 ) — - — A_(p,m 2 )J 

The zeroth order BR interaction kernel for the muonium is 

K C (E; p,q) = iR(E; p )R(E; q )V(p- q )A(p, q ), (16) 

with 

(8fxE 2 \ ^ 

(E + E p i + E p2 )(E — E p i + E p2 )(E + E pl — E p2 ) ) ^ 

We assume (0) as a definition. In the following we will verify that this choice satisfies the 
properties i) and it) given in the previous section. 

In the static limit (E — > m\ + m 2 , Ej — > rrij and p, <f — > 0), 



i^ c (E; p,q) -^-i V(p- q 



'l + 7°\ (1 7l-7 V 2)T 



2 / v 2 / 

i. e. K c reproduces the Coulomb potential times some spinorial factors. 

In order to verify that the choice (|i~6D makes solvable the BS equation (|6]) it is useful to 
express the Green function G c in terms of a new function H c 0: 

G C (E; p, q) = G (E; p) + iR{E; p )R(E; q )H C (E; p, q )G (E; p)A(p, q )G (E; q). (18) 



2 In general H c could depend on each component of the momenta p and q. The explicit calculation, 
however, will show that H c does not depend on po and go (see Eq. (p0|)). 
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Including (|T8| ) and (|T^) in (||), we have 

d 4 k R{E; k 



H c (E;p,q)=V(p-q) - i J 



(2tt) 4 (A; + E X - E kl + ie)(k° - E 2 + E k2 - it) 
x V(p-k)H c (E;k,q). 



Integrating on k° we obtain 

(27r) 3 E *_ k 2 /2 V (P~ k ) H c(&, k > V ). ( 19 ) 



where 



rf 3 A; 1 
^2tt) 3 E* - k 

(E — rrii — Tn%)(E — m\ + m 2 ){E + mi — m 2 ){E + mi + 7712) 



Equation (P~D| ) is nothing else than the Schrodinger equation for the propagator of a non 
relativistic particle in an external Coulomb field. Therefore its solution is known. A way 
to express it is by means of the Gegenbauer polynomia Cj (for the definition and some 
properties see jl3f) |I4|: 



H C {E- p, q ) = V(p- q ) - 2 ^f_J 2 , £ r— 1 _ , C) ) • f,(q )) , (20) 

where 7 = \/—2^jlE*. Substituting Eq. (|20f) in (|18|) we obtain the explicit analytic expression 
of the Green function corresponding to the kernel K c given by (|T6|) . As we will see, once G c 
is given, it is straightforward to work out the poles E^, the residua and G c n . 
From ( p0|) we have immediately that G c has poles in 



jj, fZa 
2 V~ 

Notice that up to order a 



E* = -| ( — - ) ^E = E c n = yjmi - (/iZ«/n) 2 + ^/m 2 . - (fiZa/n) 2 neN. (21) 



7-,r ^ fZaV 

E n ^m 1 +m 2 -- 1 — 1 



The poles of G c give back the mass terms plus the Bohr levels, i. e. the physically correct 
levels up to order a 2 . Moreover we have 



'mi + m 2 nil + m2 2 V n 
9 



The residuum at the pole E = E°, as defined in ([?]), is 

Y.K...(p,i) = ^ , 7 ,f, 2 7T=, ^ BMWZrt) 

^ Jmf — {jiZajnydm^ — (fiZa/n) 2 
x (p 2 + (fiZa/n) 2 ^ (q 2 + (jj,Za/n) 2 ^ r n (p, q ) 

x G (E c n ;p)A(p,q)G (E c n ;q), (22) 

where 

aA f^Za\ 5 C^_! (i(f ) • fj(q )) ' 

r n (p, q ) = 64vm — — , 2 z; 2 = LL ¥W(P )• 

V n / [p 2 + (/iZa/n) 2 ] [g 2 + (/zZa/n) 2 ] z=0 m= _, 

(23) 

(fnimip ) = Rni(p)Yim(P) are the well-known hydrogen atom wave functions flI5| . Moreover, 
we can write 

/ i + / \ « / i _ 7 o \ ( 2 ) r 

(A+(p,mi) — - — A + (g,m 1 )7°j (A_(<f,m 2 ) — - — A_(p,m 2 )7°j 

/ 1 + 7° \ / 1-7° \ 

(A+(p,mi) — - — A + (g,m 1 )7°j (A_(g,m 2 ) — - — A_(p,m 2 )7°j = 

\ - / a/3 V ~* / jS 

(A + (p>i)) ap (A + (g,m 1 )7°)^(A-(g,rn 2 )) 7 „ (^^) (A_(p, m 2 )7°) r5 

\ J pa \ / ut 

= \J2 (A + (p,m 1 )r Ss A_(p,m 2 ) 7 ) aa (A_(£ m 2 )4 8 A+(£ m 1 ) 7 °) , (24) 
where we have used the Fiertz identity: 

with the definitions: 



) = ^ E ( r LL (r Ss ) pT fife o,i, s e-fi,...,fi 

/ VT Ss 



Too = — g — 7 ' 13 = 1 — 2 — Ws ' 7 ' 

= (0,0,1) v ±1 = -^=(±l,t,0) . 

Eqs. (|23|) and (|24|) allow to identify the quantum numbers n... with the principal quantum 
number n, with the numbers S, s describing the spin of the bound state and with the numbers 
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I, m describing the angular momentum. The corresponding states RnimSs are 

z p Wmj — {j2Za/n) z Jm 2 — {jJ,Za/n) z 

i 

E pl E p2 E ql E q 2 \ 2 



x 



X 



(E pl + m 1 )(E p2 + m 2 )(E ql + m 1 )(E q2 + m 2 ) 
(p 2 + (fiZa/n) 2 ) 



x 



(p° + ^(^) - £ pl + *e)(p° - £ 2 (££) + E p2 - ie) 

(q 2 + (fiZa/n) 2 ) 

(q° + E 1 (E-) - E ql + *e)(g° - E 2 (E-) + E q2 - it) 

x (A + (p,m 1 )T Ss A-(P,m 2 h°) aS (A_(g, m 2 )4 s A+(?, ™ih°) T/3 -(25) 

These states are not degenerate (how it is possible to verify directly by calculating D n i m s s )- 
Sometimes in the literature the residua at the poles of the Green function are written as 

the functions RnimSs an d RnimSs are called the BR wave functions of the bound state: 

x / E pl E p2 \ * 

\(S p i +mi)(E p2 + m 2 )J 

(p 2 + (/iZa/n) 2 ) 

(p° + E^E*) - E pl + ie)(p° - E 2 (EZ) + E p2 - ie) 

x Vnim{p ) (A + (p,mi)r 5s A_(p,m2)7 ) a5 , (26) 

RnimSs = 1° (RnimSs) 1 1° ■ (27) 

Actually there are no reasons to introduce the bound state wave functions: from the for- 
malism developed in the previous section it is clear that all the physical quantities can be 
expressed in terms of residua. 

Finally, we obtain G c n subtracting from G c the singular part: 

G c n (E c n ,p,q) = (2,)Y(j) - g)G (i?» + Go(E»^(^;p, 9 )G (^; ?) 

+ H c (E c n ,p,q), (28) 
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where H c takes into account all the contributions non singular in E = E° which come from 
the second term in (p0|). For an explicit expression of H c in the positronium case we refer 
the reader to ||16||. 



4 Conclusions 

From the above given expressions we can recover some interesting limiting cases. 

Putting Z = 1 and mi = m 2 = m (// = m/2 and E% = E 2 = E/2) the above given 
zeroth-order solution of the BS equation reduces to the original BR solution given in || 
for positronium. The main difference with the muonium case is that for positronium also 
annihilation graphs contribute to the interaction kernel K. In the literature the singlet state 
(S = 0) is usually referred as parapositronium and the triplet state (S — 1) as orthopositro- 



nium. Some applications can be found in [12, 061 O, Qq, [L 



Taking the limit of one particle mass to infinity the case of a particle in an external 
Coulomb field is recovered. In this case (e. g. m2 — > oo ) /i = mi = m, the difference 
E2 — Tri2 is finite and the bound state energy E is given by 

E = lim (E 1 + E 2 - m 2 ). 



This case has been extensively studied in |12| ^] f° r t ne evaluation of the pure radiative 



corrections to the energy levels of hydrogenic atoms. 

As a conclusive remark, we stress that all these systems and muonium can be studied 
now in the same framework. In particular the formalism provides a powerful tool in dealing 
simultaneously with radiative, recoil and radiative-recoil corrections. Extremely interesting 
seems also to be the study of the infinit mass limit of one particle in the energy expansion 
(|T2|) evaluated on muonium states. This should eventually clarify how this limit works in a 
purely off mass-shell context. 
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